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The study of the mechanism of electroweak-symmetry 
breaking is one of the major thrusts of the experimental 
high-energy-physics program. Following the discovery of 
a Higgs-like boson at ATLAS Q and CMS [2] with a 
mass of approximately 126 GeV/c 2 and complementary 
evidence from CDF and DO [3J, the most pressing ques- 
tion is whether this state is in fact the Higgs boson of the 
standard model (SM), part of an extended Higgs sector 
(such as that of the minimal supersymmctric standard 
model, MSSM 01), a composite Higgs or a completely 
different particle with Higgs-like couplings (such as a ra- 
dion in warped extra dimensions @ or a dilaton Q). 

We search for particles in an extension to the standard 
model that includes a light neutral Higgs boson, h°, with 
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mass m^o = 126 GeV/c 2 . Rather than assuming a par- 
ticular theoretical framework (such as the MSSM), we 
follow a phenomenological approach, using a general two- 
Higgs doublet model as a convenient simplified model Q , 
which contains a heavy charged Higgs boson and a 
heavier neutral state H°. In this approach, the search for 
a number of specific final states that have the strongest 
couplings to Higgs particles is motivated 0, [Io| • The fi- 
nal state of a W-boson pair (WW) is enhanced by WW 
scattering in models where the Higgs sector is strongly 
coupled [ll|. This signal has been the subject of much 
detailed investigation [l2|- The phenomenology of reso- 
nant production of the final states Zh° [IH and W + W~ Z 
fl4| has also been investigated. 

In this letter, we focus on the final state W + W~bb [lij ]. 
which can have a large production rate from the pro- 
cess gg H° followed by H° i/ ± W T with ff± 
W ± h° ->■ WHb. The W+W~bb final state is also the 
final state of top-quark pair production, and has been 
extensively studied. However, no search for Higgs-boson 
cascades as described here has been reported previously, 
though searches have been performed for charged Higgs 
bosons in top-quark decays t — > H ± b IB- 

We analyze a data sample corresponding to an inte- 
grated luminosity of 8.7±0.5 fb _1 recorded by the CDF 
II detector [lij], a general purpose detector designed to 
study pp collisions at = 1.96 TeV in the Fermilab 
Tevatron collider. The CDF tracking system consists of 
a silicon microstrip tracker and a drift chamber that are 
immersed in a 1.4 T axial magnetic field [lo[. Projective- 
tower-geometry electromagnetic and hadronic calorime- 
ters surrounding the tracking system measure particle en- 
ergies, with muon detection provided by additional drift 
chambers located outside the calorimeters. 

The signature of H° ->■ W =F H ± W'W+h 
W~~W + bb is a charged lepton (e or n), missing trans- 
verse momentum, two jets arising from b quarks, and two 
additional jets from a W-boson hadronic decay. Events 
are selected online (triggered) by the requirement of an 
electron (e) or muon (fj,) candidate [2l| with transverse 
momentum pt 22j greater than 18 GeV/c. After trig- 
ger selection, events are retained if the electron or muon 
candidate has a pseudorapidity \rj\ < 1.1 [22J, px > 20 
GeV/c, and satisfies the standard CDF identification and 
isolation requirements 21] . We reconstruct jets in the 
calorimeter using the jetclu [23} algorithm with a clus- 
tering radius of 0.4 in 77 — <j> space. The jets are cal- 



ibrated using the techniques outlined in Ref. [24|. At 
least four jets are required, each with transverse energy 
Et > 15 GeV and |?7| < 2.4. Missing transverse momen- 
tum [25| is reconstructed using calorimeter and muon in- 
formation 21] ; in the W + W~bb experimental signature, 
the missing transverse momentum is mostly due to the 
neutrino from the leptonically-decaying W boson. We 
require Et > 20 GeV/c. Since such a signal would yield 
two jets originating from b quarks, we require (with min- 



imal loss of efficiency) evidence of decay of a b hadron 
in at least one jet. This requirement, called 6-tagging, 
makes use of the SECVTX algorithm, which identifies jets 
from b quarks via their secondary vertices [26j |. 

We model the production of H° bosons with mjo = 
325-1100 GeV/c 2 and subsequent decays H° -> W^H ± 
with m H ± = 225-600 GeV/c 2 and decays H ± -> W ± h° 



23. Ad- 



with m h o = 126 GeV/c 2 , all with madgraph 
ditional radiation, hadronization, and showering are de- 
scribed by pythia [28]. The detector response for all 
simulated samples is modeled by the GEANT-based CDF 
II detector simulation [29I ]. 

The dominant SM background to this signature is top- 
quark pair production. We model this background using 
pythia with a top-quark mass m t = 172.5 GeV/c 2 [3p| | . 
We normalize the ti background to the theoretical cal- 
culation at next-to-next-to-leading order (NNLO) in the 
strong interaction coupling constant, a s (3lj. In addi- 
tion, events generated by a next-to-leading order pro- 
gram, MC@NLO [32I are used in estimating an uncertainty 
in modeling the radiation of an additional jet. 

The second largest SM background process is the as- 
sociated production of a W boson and jets. Samples of 
W-boson+jets events with light- and heavy-flavor (6, c) 
quark jets are generated using ALPGEN [331 ] . and inter- 
faced with a parton-shower model from PYTHIA. The W- 
boson+jets samples are normalized to the measured W- 
boson-production cross section, with an additional mul- 
tiplicative factor for the relative contribution of heavy- 
and light-flavor jets, following Ref. [2f|. 

Backgrounds due to production of a Z boson with 
additional jets, where the second lepton from the Z- 
boson decay is not reconstructed, are small compared to 
the W-boson background and are modeled using events 
generated with ALPGEN interfaced to the parton-shower 
model from pythia. The multi-jet background, in which 
a jet is misreconstructed as a lepton, is modeled using 
events triggered on jets and normalized to a background- 
dominated region at low missing transverse momentum 
where the multi-jet background is large. 

The SM backgrounds due to production of single top 
quarks and pairs of vector bosons are modeled using 
MADGRAPH interfaced with pythia parton-shower mod- 
els and pythia, respectively, and normalized to next-to- 
leading-order cross sections [34l. l35|. 

The Higgs-boson candidate mass reconstruction begins 
with identification of the leptonically-decaying W boson, 
assuming the missing transverse momentum is due to the 
resulting neutrino. Of the multiple solutions for the neu- 
trino pseudorapidity, we use the smallest value that yields 
the reconstructed W mass closest to the known value. 
The hadronically-decaying W boson is identified as the 
pair of jets that yield the reconstructed dijet mass clos- 
est to the known W mass, excluding jets with a 6-tag. If 
fewer than two jets without 6-tags are present, the same 
procedure is used but modified to include the 6-tagged 
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FIG. 1: Distribution of reconstructed Higgs-boson masses in 
simulated events. Top: m h a = 126 GeV/c 2 reconstructed as 
mtb, center: m H ± as mwbb, and bottom: m H o as rriwwbb- 



by theoretical uncertainties on the cross sections of the 
background processes. To probe the description of ad- 
ditional jets, we compare our nominal ti model to one 
generated by mc@nlo and take the full difference as a 
systematic uncertainty. We also consider systematic un- 
certainties associated with the description of initial- and 
final-state radiation }36J , uncertainties in the efficiency of 
reconstructing leptons and identifying 6-quark jets, and 
uncertainties in the contribution from multiple interac- 
tions. In addition, we consider a variation of the Q 2 scale 
of W-boson+jet events in ALGPEN. In each case, we treat 
the unknown underlying quantity as a nuisance parame- 
ter. Except in the case of the normalization uncertainty, 
which affects only the overall rates, for each source of un- 
certainty we measure the distortion of the m b i spectrum 
for positive and negative fluctuations of the underlying 
quantity. Table U lists the contributions of each of these 
sources of systematic uncertainty to the yields. 



TABLE I: Contributions to the systematic uncertainty on the 
expected numbers of events for the two main background 
processes, the total background yield, and an example 500 
GeV/c 2 Higgs-boson signal with an assumed total cross sec- 
tion of 1 pb. 



Process 


ti 


W-boson+jets 


Total bg. 


Higgs 


Predicted yield 


229 


43 


294 


341 


Jet energy scale 


23% 




17% 


12% 


Radiation 


3% 




2% 


8% 


Q 2 scale 




18% 


3% 




Mult, interactions 


1% 


6% 


2% 




ti generator 


5% 




4% 




Normalization 


10% 


30% 


16% 




Total syst. uncert. 


26% 


35% 


24% 


15% 



jets. The light hP is reconstructed from the remaining 
6-tagged jets. If fewer than two 6-tagged jets remain, 
the jet or jets with largest transverse momentum not 
associated with the hadronic VF-boson decay are used 
instead, without significant loss of mass resolution. Fig- 
ure [1] shows distributions of the reconstructed mass for 
several choices of Higgs masses. 

We enhance the signal-to-background ratio through re- 
quirements on the mass of the W + W~bb and W^bb sys- 
tems, and search for an excess of events above expec- 
tations from backgrounds in event distributions versus 
the mass of the bb system (h° —> bb). Backgrounds have 
broad, smoothly decreasing distributions while a signal 
would be reconstructed near the Higgs-boson mass. 

We consider several sources of systematic uncertainty 
on the predicted background rates and distributions, as 
well as on the expectations for a signal. Each systematic 
uncertainty affects the expected sensitivity to a signal, 
expressed as an expected cross-section upper limit in the 
no-signal assumption. The dominant systematic uncer- 
tainty is the jet-energy-scale uncertainty 24|, followed 



We validate our modeling of the SM backgrounds in 
four background-dominated control regions. Each con- 
trol region preserves the one lepton and at least four 
jet requirements with additional requirements per region. 
Events in the first region are used to study the W + W~bb 
and W ± bb mass reconstruction, requiring at least one b- 
tagged jet and bb mass smaller than 100 GeV/c 2 . The 
second region probes bb and W + W~bb mass reconstruc- 
tion, requiring at least one 6-tagged jet and W ± bb mass 
smaller than 250 GeV/c 2 . The third region tests the 
modeling of W ± bb and bb mass reconstruction, requiring 
at least one 6-tagged jet and W + W~bb mass less than 
450 GeV/c 2 . The fourth region tests the modeling of the 
VF-boson-plus-jets background, requiring exactly zero b- 
tagged jets and W + W~bb mass greater than 450 GeV/c 2 . 
Assuming an H° production cross section of 250 fb, each 
control region is expected to have negligible signal con- 
tamination, with the exception of the zero 6-tag region 
which would include signal events at approximately 10% 
of the sample. For two of the control regions, Fig.[2]shows 
the reconstructed bb mass distributions which, along with 
other similar distributions, indicate that the background 
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FIG. 2: Distribution of events versus reconstructed bb invari- 
ant mass (m b i) for observed data and expected backgrounds 
in two control regions. Top, control region consisting events 
with at least four jets, exactly zero 6-tags and mwwbb < 450 
GeV/c 2 . Bottom, control region consisting of events with 
at least four jets and mwbb < 250 GeV/c 2 . The lower panels 
give the relative difference between the observed and expected 
distributions; the hatched areas show the combined statisti- 
cal and systematic uncertainties of the expected background. 
The small dip near 80 GeV/c 2 is mainly due to the VK-boson 
mass reconstruction. 
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FIG. 3: Distribution of events versus reconstructed bb invari- 
ant mass (m b i), for observed data and expected backgrounds 
in the signal region. A signal hypothesis is shown, assum- 
ing a total cross section of 250 fb, m H o — 500 GeV/c 2 , and 
m H ± = 300 GeV/c 2 . See Fig[2]for descriptions of lower panel 
and hatching. 



blcEO 




FIG. 4: Upper limits at 95% C.L. on the cross section times 
branching fraction as a function of the Higgs-boson masses 
m H ± and m H o; m h o is fixed to 126 GeV/c 2 in each case. 
Diamonds show the grid of probed masses; the intermediate 
values are interpolated. 



mass distributions are well modeled within systematic 
uncertainties. 

Figure [3] shows the observed distribution of events in a 
representative signal region compared to possible signals 
and estimated backgrounds. At each Higgs-boson mass 
hypothesis, we fit the most likely value of the Higgs-boson 
cross section by performing a maximum-likelihood fit in 
the binned m b i distribution, allowing for systematic and 
statistical fluctuations via template morphing [37|. No 
evidence is found for the presence of Higgs-boson cas- 
cade decays in WWbb events. We set upper limits on 
Higgs production at 95% confidence level using the CLs 
method [38| . without profiling the systematic uncertain- 
ties. The observed limits are consistent with expectation 
for the background-only hypothesis. See Fig. @] and Ta- 



in conclusion, we report on the first search for multi- 
ple Higgs bosons in cascade decays. For each accepted 
event, we reconstruct the lightest neutral Higgs-boson 
mass (jn b i), and find the CDF data to be consistent with 
standard model background predictions. We calculate 
95% C.L. upper limits on the cross section of such Higgs- 
boson production, assuming 100% branching ratio of H° 
to W ± HT and H* to W ± h°, from 1.3 pb to 0.015 pb for 
masses ranging from (m H o = 325, m H ± = 225) GeV/c 2 
to (m H a = 1100, rn#± = 600) GeV/c 2 respectively, and 
interpret the limits in terms of a simplified two-Higgs- 
doublet model. While the limits cited here do not ex- 
clude any region in the m#o — m H ± -plane in the simpli- 
fied model used, there are the first such limits available. 
The larger center-of-mass energy and integrated luminos- 
ity of data collected by the LHC experiments are likely to 
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TABLE II: Signal region definitions and expected and ob- 
served 95% C.L. upper limits on the production cross section 
times branching fraction for each Higgs-boson mass hypothe- 
sis. Theoretical predictions are also shown [39h41|. 



(m H o,m H ±) 






Exp (Obs) 


Theory 


(GoV/c 2 ) 


(GoV/c 2 ) 


(GoV/c 2 ) 


Limit (fb) 


(fb) 


325, 225 


> 175 


> 275 


1100 (1300) 


34 


400, 300 


> 225 


> 325 


960 (1100) 


18 


425, 225 


> 200 


> 375 


900 (960) 


13 


500, 300 


> 200 


> 450 


470 (590) 


3.9 


500, 400 


> 350 


> 450 


510 (700) 


3.9 


525, 225 


> 100 


> 500 


420 (460) 


2.5 


DUU, OUU 


~> ZUU 


OOU 


900 flSO^ 


7fi 


600, 400 


> 350 


> 550 


210 (250) 


0.76 


700, 400 


> 325 


> 650 


90 (100) 


0.15 


700, 600 


> 450 


> 650 


10 (96) 


0.15 


725, 225 


> 425 


> 700 


90 (120) 


0.10 


800, 300 


> 275 


> 750 


50 (51) 


3 x 10~ 2 


800, 600 


> 475 


> 725 


43 (46) 


3 x 10~ 2 


900, 400 


> 450 


> 775 


28 (36) 


6 x 10 -3 


900, 600 


> 475 


> 800 


24 (29) 


6 x 10" 3 


1100, 600 


> 475 


> 975 


13 (15) 


2 x 10" 4 



have the sensitivity to discover or exclude such models. 
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